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Abstract 

Background: Malaria caused by Plasmodium vivax is still a public health problem in the Republic of Korea (ROK), 
particularly regarding the recent re-emergence of this malarial species near the demilitarized zone in northwestern 
Paju City, Gyeonggi-do Province. Currently, at least 4 species {An. kleini, An. pullus, An. belenrae and An. lesteri) of the 
Hyrcanus Group are reported as possible natural vectors of vivax malaria in the ROK, and An. sinensis, which is the 
most dominant species, has long been incriminated as an important natural vector of this P. vivax. However, An. 
sinensis was ranked recently as a low potential vector. According to the discovery of natural hybrids between An. 
sinensis (a low potential vector for P. vivax) and An. kleini (a high potential vector for P. vivax) in Paju City, intensive 
investigation of this phenomenon is warranted under laboratory conditions. 

Methods: Mosquitoes were collected during 2010-2012 from Paju City, ROK. Hybridization experiments used 
iso-female line colonies of these anophelines together with DNA analysis of ribosomal DNA [second internal 
transcribed spacer (ITS2)] and mitochondrial DNA [cytochrome c oxidase subunit I (COI)] of the parental colonies, 
F r hybrids and repeated backcross progenies were performed intensively by using a PCR-based assay and 
pyrosequencing technology. 

Results: The results from hybridization experiments and molecular investigations revealed that the mitochondrial 
COI gene was introgressed from An. sinensis into An. kleini. The An. sinensis progenies obtained from consecutive 
repeated backcrosses in both directions, i.e., F 2 -n progeny [{An. sinensis xAn. kleini) xAn. sinensis] and F 3 - 5 progeny 
[{An. kleini x An. sinensis) x An. kleini] provided good supportive evidence. 

Conclusions: This study revealed introgression of the mitochondrial COI gene between An. sinensis and An. kleini 
through consecutive repeated backcrosses under laboratory conditions. This new body of knowledge will be 
emphasized in reliable promising strategies in order to replace the population of An. kleini as a high potential 
vector for P. vivax, with that of a low potential vector, An. sinensis, through the mechanism of gene introgression in 
nature. 

Keywords: Anopheles sinensis, An. kleini, Hybridization experiment, Second internal transcribed spacer, Cytochrome c 
oxidase subunit I, Introgression 
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Background 

Up until now, at least 26 species members of the Anopheles 
hyrcanus group have been reported, and their distribution 
has extended widely from Europe to East and Southeast 
Asia, including some of the off-lying islands of the Indian 
and Pacific Oceans [1]. Some species of the Hyrcanus 
Group are accepted as important vectors in transmitting 
human diseases, e.g., malaria {Plasmodium vivax) [2-14], 
filariae {Wuchereria bancrofti and Brugia malayi) [15,16], 
and Japanese encephalitis virus [17,18], particularly in 
the Oriental and contiguous parts of eastern Palaearctic 
regions. 

At least 6 species {Anopheles sinensis, An. lesteri, An. 
pullus, An. sineroides, An. belenrae and An. kleini) of the 
Hyrcanus Group are found indigenously in the Republic 
of Korea (ROK). Among these, An. sinensis has long been 
incriminated as the most dominant and important natural 
vector of P. vivax, especially due to the recent re- 
emergence of vivax malaria near the demilitarized zone in 
northwestern Paju City, Gyeonggi-do Province [3,19-22]. 
However, the low concentration of circumsporozoite (CS) 
antigen obtained from wild-caught females [6], and very 
low sporozoite rates recovered from laboratory suscepti- 
bility tests [7,9], have brought about the cryptic status of 
An. sinensis as a natural vector of vivax malaria transmis- 
sion in the ROK. Consequently, the implication of other 
An. hyrcanus species, i.e., An. kleini, An. pullus, An. belen- 
rae and An. lesteri as possible natural vectors of vivax mal- 
aria in the ROK has been proposed extensively [8,9], even 
though the latter species is thought to have a small popu- 
lation [7]. Remarkably, An. sinensis strain from China 
has been incriminated recently as an efficient vector of 
P. vivax [11]. 

The discovery of natural hybrids has been reported 
from some important malaria vectors, for instance, be- 
tween An. gambiae and An. arabiensis [23], An. scanloni 
{= dims C) and An. baimaii {= dims D) [24], and An. 
minimus and An. harrisoni [25,26]. Regarding An. kleini 
(a high potential vector for P. vivax) and An. sinensis 
(a low potential vector for P. vivax) [7,9], a single gravid, 
natural hybrid female was discovered between these 2 
anopheline species from a total of 658 wild-caught females 
in Paju City, ROK; an endemic zone of vivax malaria from 
2004 to 2008. These wild-caught females comprised 360 
An. sinensis, 258 An. pullus, 20 An. belenrae, 15 An. kleini, 
3 An. sineroides, 1 An. lesteri and 1 hybrid female between 
An. kleini and An. sinensis [27]. The results of self-crossing 
between F x -progenies derived from one egg-batch of an 
egg-laid gravid female, and molecular investigations, 
revealed possible natural backcrossing (introgression) 
between a hybrid female and male An. sinensis. There is 
no clear basic information on the introgression between 
An. sinensis and An. kleini, as proposed by [27]. Thus, 
systematic investigations into the role of introgressive 



hybridization between these two anopheline species were 
performed in this study. Accordingly, attempts were made 
to establish iso-female line colonies of An. sinensis and 
An. kleini, perform crossing experiments (reciprocal and 
back crosses), investigations of reproductive systems of 
hybrid and backcross progenies, and compare DNA ana- 
lysis of ribosomal DNA (ITS2) and mitochondrial DNA 
(COI) of the parental, F^hybrids with that of repeated 
backcross progenies by using a PCR-based assay and pyro- 
sequencing technology. 

Methods 

Establishment of iso-female lines 

Mosquitoes were collected during 2010-2012 from Paju 
City, Republic of Korea (ROK), by placing a light trap in 
cowsheds. Then, wild-caught females were transported 
for colonization to the insectary of the Department of 
Parasitology, Faculty of Medicine, Chiang Mai University, 
Chiang Mai, Thailand. Four iso-female lines of both An. 
sinensis and An. kleini were established successfully using 
the methods of [28]. An Fi-progeny of each iso-female line 
was used for species identification following the keys of 
[29] as well as a molecular assay [30]. Then, one iso- 
female line of each species, with molecular identification 
of both nuclear (ITS2) and mitochondrial (COI) genes, 
were well matched with those in the GenBank nucleotide 
sequence database, and selected, i.e., An. sinensis F 0 -l 
(SF0-1) and An. kleini F 0 -l (KF0-1). These iso-female lines 
have been maintained in colonies in the laboratory at 
Chiang Mai University for more than 10 consecutive gen- 
erations, and used for hybridization experiments and com- 
parative DNA sequence analyses. 

Hybridization experiments 

Hybridization experiments (reciprocal and back crosses, 
and repeated backcross progenies) between An. sinensis 
and An. kleini were performed by using virgin females 
and males and following the techniques previously re- 
ported by [31]. Post-mating reproductive isolation was 
recorded using the criteria of low viability (hatchability, 
survival, pupation, and emergence), adult sex distortion 
and abnormal morphology of the reproductive system. 

PCR identification, dideoxy sequencing and 
phylogenetic analysis 

DNA was extracted individually from 60 mosquitoes using 
the RED Extract-N-Amp™ Tissue PCR kit ( Sigma- Aldrich, 
Spruce Street, SL) as shown in Table 1. Primers for the 
amplification of ITS2 and COI regions followed a previous 
report by [30]. The ITS2 region of the rDNA was ampli- 
fied using the ITS2 Forward and ANO 28S-20 primers 
[30,32]. The mitochondrial COI gene was amplified using 
the LCO1490 (5 -GGT CAA CAA ATC ATA AAG ATA 
TTG G-3 ) and HC02198 (5 -TAA ACT TCA GGG 
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Table 1 Species identification of An. sinensis and An. 
GenBank accession numbers 



kleini samples based on ITS2 and COI sequences, and their 



Samples (Female x Male) 



Code of samples Identified species name (GenBank accession number) 

ITS2 COI 



Reference 



Parental 

An. sinensis F 0 -1 

An. sinensis F 0 -2 

An. sinensis F 0 -3 

An. sinensis F 0 -4 

An. kleini F 0 -1 

An. kleini F 0 -2 

An. kleini F 0 -3 

An. kleini F 0 -4 

Reciprocal crosses 

An. sinensis x An. kleini F r 1 

An. sinensis x An. kleini Ft -2 

An. sinensis x An. kleini Ft -3 

An. kleini x An. sinensis F r 1 

/An. /c/e/n/ x An. sinensis Ft -2 

/\n. /c/e/n/' x An. sinensis F r 3 

Back crosses 

(/\n. sinensis x /An. /c/e/n/) x /An. sinensis 



SFO-1 
SFO-2 
SFO-3 
SFO-4 
KFO-1 
KFO-2 
KFO-3 
KFO-4 

SKF1-1 
SKF1-2 
SKF1-3 
KSF1-1 
KSF1-2 
KSF1-3 



An. sinensis (KC797396) 
An. sinensis (KC797397) 
An. sinensis (KC797398) 
An. sinensis (KC797399) 
An. kleini (KC797431) 
An. kleini (KC797432) 
An. kleini (KC797433) 
An. kleini (KC797434) 

Mixed 
Mixed 
Mixed 
Mixed 
Mixed 
Mixed 



An. sinensis (KC797435) 
An. sinensis (KC797436) 
An. sinensis (KC797437) 
An. sinensis (KC797438) 
An. kleini (KC797439) 
An. sinensis (KC797440) 
An. sinensis (KC797441) 
An. sinensis (KC797442) 

An. sinensis (KC797446) 
An. sinensis (KC797447) 
An. sinensis (KC797448) 
An. sinensis (KC797443) 
An. kleini (KC797444) 
An. sinensis (KC797445) 



This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 

This study 
This study 
This study 
This study 
This study 
This study 



hybridFq 


x sinensis-] 


= > (hybridF 2 


-1) 


SKSF2-1 




Mixed 


An. 


sinensis 


(KC797449) 


This study 


hybridFq 


x sinensis-2 


= > (hybridF 2 


-2) 


SKSF2-2 




Mixed 


An. 


sinensis 


(KC797450) 


This study 


hybridFq 


x sinensis-?) 


= > (hybridF 2 


-3) 


SKSF2-3 




Mixed 


An. 


sinensis 


(KC797451) 


This study 


hybridF 2 


x sinensis-] 


= > (hybridF 3 


-D 


SKSF3-1 


An. 


sinensis (KC797400) 


An. 


sinensis 


(KC797452) 


This study 


hybridF 2 


x sinensis-2 


= > (hybridF 3 


-2) 


SKSF3-2 


An. 


sinensis (KC797401) 


An. 


sinensis 


(KC797453) 


This study 


hybridF 2 


x sinensis-? 


= > (hybridF 3 


-3) 


SKSF3-3 


An. 


sinensis (KC797402) 


An. 


sinensis 


(KC797454) 


This study 


hybridF 3 


x sinensis-] 


= > (hybridF 4 


-D 


SKSF4-1 


An. 


sinensis (KC797403) 


An. 


sinensis 


(KC797455) 


This study 


hybridF 3 


x sinensis-2 


= > (hybridF 4 


-2) 


SKSF4-2 


An. 


sinensis (KC797404) 


An. 


sinensis 


(KC797456) 


This study 


hybridF 3 


x sinensis-? 


= > (hybridF 4 


-3) 


SKSF4-3 


An. 


sinensis (KC797405) 


An. 


sinensis 


(KC797457) 


This study 


hybridF 4 


x sinensis-] 


= > (hybridF 5 


-D 


SKSF5-1 


An. sinensis (KC7974006) 


An. 


sinensis 


(KC797458) 


This study 


hybridF 4 


x sinensis-2 


= > (hybridF 5 


-2) 


SKSF5-2 


An. 


sinensis (KC797407) 


An. 


sinensis 


(KC797459) 


This study 


hybridF 4 


x sinensis-? 


= > (hybridF 5 


-3) 


SKSF5-3 


An. 


sinensis (KC797408) 


An. 


sinensis 


(KC797460) 


This study 


hybridF 5 


x sinensis-] 


= > (hybridF 6 


-D 


SKSF6-1 


An. 


sinensis (KC797409) 


An. 


sinensis 


(KC797461) 


This study 


hybridF 5 


x sinensis-2 


= > (hybridF 6 


-2) 


SKSF6-2 


An. 


sinensis (KC797410) 


An. 


sinensis 


(KC797462) 


This study 


hybridF 5 


x sinensis-? 


= > (hybridF 6 


-3) 


SKSF6-3 


An. 


sinensis (KC797411) 


An. 


sinensis 


(KC797463) 


This study 


hybridF 6 


x sinensis-] 


= > (hybridF 7 


-D 


SKSF7-1 


An. 


sinensis (KC797412) 


An. 


sinensis 


(KC797464) 


This study 


hybridF 6 


x sinensis-2 


= > (hybridF 7 


-2) 


SKSF7-2 


An. 


sinensis (KC797413) 


An. 


sinensis 


(KC797465) 


This study 


hybridF 6 


x sinensis-? 


= > (hybridF 7 


-3) 


SKSF7-3 


An. 


sinensis (KC797414) 


An. 


sinensis 


(KC797466) 


This study 


hybridF 7 


x sinensis-] 


= > (hybridF 8 


-D 


SKSF8-1 


An. 


sinensis (KC797415) 


An. 


sinensis 


(KC797467) 


This study 


hybridF 7 


x sinensis-2 


= > (hybridF 8 


-2) 


SKSF8-2 


An. 


sinensis (KC797416) 


An. 


sinensis 


(KC797468) 


This study 


hybridF 7 


x sinensis-? 


= > (hybridF 8 


-3) 


SKSF8-3 


An. 


sinensis (KC797417) 


An. 


sinensis 


(KC797469) 


This study 


hybridF 8 


x sinensis-] 


= > (hybridF 9 


-D 


SKSF9-1 


An. 


sinensis (KC797418) 


An. 


sinensis 


(KC797470) 


This study 


hybridF 8 


x sinensis-2 


= > (hybridF 9 


-2) 


SKSF9-2 


An. 


sinensis (KC797419) 


An. 


sinensis 


(KC797471) 


This study 


hybridF 8 


x sinensis-? 


= > (hybridFg 


-3) 


SKSF9-3 


An. 


sinensis (KC797420) 


An. 


sinensis 


(KC797472) 


This study 
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Table 1 Species identification of An. sinensis and An. kleini samples based on ITS2 and COI sequences, and their 
GenBank accession numbers (Continued) 



hybridF 9 x sinensis-] = > (hybridF 10 -1) 


SKSF10-1 


An. sinensis (KC797421) 


An. sinensis (KC797473) 


This study 


hybrid F 9 x sinensis-2 = > (hybrid F ] 0 -2) 


SKSF10-2 


An. sinensis (KC797422) 


An. sinensis (KC797474) 


This study 


hybrid F 9 x sinensis-3 = > (hybrid Ft 0 -3) 


SKSF10-3 


An. sinensis (KC797423) 


An. sinensis (KC797475) 


This study 


hybrid Ft o x sinensis-] = > (hybridF n -1) 


SKSF1 1-1 


An. sinensis (KC797424) 


An. sinensis (KC797476) 


This study 


hybrid Ft o x sinensis-2 = > (hybrid F ] ] -2) 


SKSF11-2 


An. sinensis (KC797425) 


An. sinensis (KC797477) 


This study 


hybrid Ft o x sinensis-3 = > (hybrid Ft } -3) 


SKSF11-3 


An. sinensis (KC797426) 


An. sinensis (KC797478) 


This study 


hybridFn x hybridF n -1 


HF11-1 


An. sinensis (KC797427) 


An. sinensis (KC797479) 


This study 


hybridF n x hybridF n -2 


HF11-2 


An. sinensis (KC797428) 


An. sinensis (KC797480) 


This study 


hybridF n x hybridF n -3 


HF11-3 


An. sinensis (KC797429) 


An. sinensis (KC797481) 


This study 


hybridFn x hybridF ir 4 


HF11-4 


An. sinensis (KC797430) 


An. sinensis (KC797482) 


This study 


Back crosses 










(An. kleini x An. sinensis) x An. kleini 










hybridF] x kleini-] = > (hybridF 2 -1) 


KSKF2-1 


An. kleini (KC890843) 


An. kleini (KC797483) 


This study 


hybrid F] x kleini-2 = > (hybrid F 2 -2) 


KSKF2-2 


An. kleini (KC890844) 


An. kleini (KC797484) 


This study 


hybrid F] x kleini-3 = > (hybrid F 2 -3) 


KSKF2-3 


An. kleini (KC890845) 


An. kleini (KC797485) 


This study 


hybridF 2 x kleini-] = > (hybridF 3 -1) 


KSKF3-1 


An. kleini (KC890846) 


An. sinensis (KC797486) 


This study 


hybrid F 2 x kleini-2 = > (hybrid F 3 -2) 


KSKF3-2 


An. kleini (KC890847) 


An. sinensis (KC797487) 


This study 


hybrid F 2 x kleini-3 = > (hybrid F 3 -3) 


KSKF3-3 


An. kleini (KC890848) 


An. sinensis (KC797488) 


This study 


hybridF 3 x kleini-] = > (hybridF 4 -1) 


KSKF4-1 


An. kleini (KC890849) 


An. sinensis (KC797489) 


This study 


hybrid F 3 x kleini-2 = > (hybrid F 4 -2) 


KSKF4-2 


An. kleini (KC890850) 


An. sinensis (KC797490) 


This study 


hybrid F 3 x kleini-3 = > (hybrid F 4 -3) 


KSKF4-3 


An. kleini (KC890851) 


An. sinensis (KC797491) 


This study 


hybridF 4 x kleini-] = > (hybridF 5 -1) 


KSKF5-1 


An. kleini (KC890852) 


An. sinensis (KC797492) 


This study 


hybrid F 4 x kleini-2 = > (hybrid F 5 -2) 


KSKF5-2 


An. kleini (KC890853) 


An. sinensis (KC797493) 


This study 


hybrid F 4 x kleini-3 = > (hybrid F 5 -3) 


KSKF5-3 


An. kleini (KC890854) 


An. sinensis (KC797494) 


This study 


An. sinensis 




An. sinensis (GU384700) 




[30] 








An. sinensis {AY 444351) 


[39] 


An. kleini 




An. kleini (GU 3847 13) 




[30] 








An. kleini (GQ265917) 


[27] 



TGA CCA AAA AAT CA-3') primers of [33]. PCR reac- 
tion was carried out in a total volume of 25 \A containing 
10 pM of each primer; and 2.5 \A of 10X buffer containing 
50 mM KC1, 10 mM Tris-HCI, 0.1% Triton X 100 supple- 
mented with 1.5 mM MgCl 2 (Promega, USA), 200 uM of 
each dNTP (GeneCraft, Germany), 0.5 \A of Taq DNA 
polymerase (Promega, USA) and 10-100 pg of genomic 
DNA. The amplification profile comprised initial denatur- 
ation at 94°C for 3 min, with 30 cycles at 94°C for 30 sec, 
55°C for 30 sec, and 72°C for 2 min, and a final extension 
at 72°C for 7 min. The PCR products were separated by 
electrophoresis on a 1.5% agarose gel stained with eth- 
idium bromide. Finally, the purified PCR products were 
subjected to sequencing in an ABI PRISM 3700 DNA 
Analyzer (Applied Biosystems, Foster City, CA) using a 
Dye Terminator Cycle Sequencing Ready Reaction Kit 
(Applied Biosystems). The sequence data obtained were 
deposited in the GenBank nucleotide sequence database 



(Table 1). The newly determined ITS2 and COI sequences 
were also compared with those available in GenBank, 
using the Basic Local Alignment Search Tool (BLAST) 
available at http://blast.ncbi.nlm.nih.gov/Blast.cgi. The 
DNA sequence data were edited manually in BioEdit ver- 
sion 7.0.5.3 [34] and aligned using CLUSTAL W [35]. 
Constructions of neighbor- joining trees using distance 
[36], and the bootstrap test with 1,000 replications, were 
performed with the MEGA version 4.0 program based on 
COI sequences [37]. The COI sequences of An. peditae- 
niatus (GenBank accession number AB539069) [38] and 
An. pullus (GenBank accession number AB444348) [39] 
were included in phylogenetic analysis. 

PCR and pyrosequencing for the detection of An. sinensis 
and An. kleini sequences 

DNA was extracted from each adult female mosquito of 
An. sinensis and An. kleini using a NucleoSpin tissue kit 



Choochote et al. Parasites & Vectors 2014, 7:36 
http://www.parasitesandvectors.eom/content/7/1/36 



Page 5 of 13 



(Macherey-Nagel GmbH and Co., Duren, Germany). 
The procedure for the pyrosequencing assay followed 
that previously described by [40]. The forward primer 
(Anop_COI_F: 5 -GAG CCC CTG ATA TAG CTT 
TTC CT-3'), and biotinylated reverse primers (Anop_- 
COI_Rb: 5'- Biotin-CCA GAT GAA AGT GGG GGA 
TAA -3'), were designed to amplify a 142-bp fragment 
of COI, and a primer Anop_COI_S (5' -ATA AGT TTT 
TGA ATA TTA CC -3') for pyrosequencing. Positive- 
control plasmids of each species were constructed by 
amplification of the 142-bp PCR products using Anop_- 
COI_F and Anop_COI_R primers. They were ligated 
and transformed into a pGEM -T Easy vector (Promega, 
WI) and an Escherichia coli JM109, respectively. The re- 
combinant plasmids were sequenced bidirectionally in 
order to confirm the correction of data. The 142 bp was 
amplified from genomic DNA using the Anop_COI_F 
and Anop_COI_Rb primers. The reactions of PCR amp- 
lification were performed in a total volume of 25 [A con- 
taining IX PCR buffer (Invitrogen, Carlsbad, CA) with 
0.2 mM of each dNTP, 2 mM MgS04, 0.4 (iM of each 
primer, 0.625 U of Platinum Taq DNA polymerase high 
fidelity (Invitrogen, Carlsbad, CA) and 2 [A of the DNA 
sample. The PCR assay was conducted using a Gen- 
eAmp PCR system 9700 thermal cycler (Applied Biosys- 
tems, Singapore). The thermocycling program included: 
initial denaturation at 94°C for 5 min, with 40 cycles at 
94°C for 30 sec, 57°C for 30 sec, and 72°C for 30 sec, 
and a final extension at 72°C for 7 min. PCR amplicons 
were detected by electrophoresis on 1.5% agarose gels. 
For pyrosequencing assays, 20 \A of the biotinylated PCR 
product of each sample was immobilized in the binding 
buffer with Streptavidin Sepharose™ beads (GE Health- 
care Biosciences AB, Uppsala, Sweden). The beads to- 
gether with the DNA were aspirated to a 96 format filter 
tool and passed through 70% ethanol and 0.2 M NaOH, 
and then washed with 10 mM Tris-acetate (pH 7.6) 
using a PyroMark™ Vacuum Prep Workstation (Biotage, 
Uppsala, Sweden). The beads were released subsequently 
into a PSQ™ 96 plate low (Biotage) containing 0.4 \iM 
Anop_COI _S sequencing primer in the annealing buf- 
fer. The samples were heated to 80°C for 2 min before 
performing pyrosequencing reactions using PyroMark™ 
Gold Q96 SQA reagents and subjecting them to the 
PyroMark™ Q96 ID instrument (Biotage). Positive and 
negative controls were included in each pyrosequencing 
assay. Finally, the PyroMark™ Q96 ID software version 
1.0 was used to analyze the results and generate a 
pyrogram. 

Results 

Hybridization experiments 

The hatchability, pupation, emergence and adult sex- 
ratio of parental, reciprocal and back crosses, repeated 



backcrosses and hybrid crosses between An. sinensis and 
An. kleini were 88.06-92.57%, 77.96-92.97%, 94.10-96.98% 
and 0.81-0.89; 82.94-85.07%, 98.00-100%, 98.87-100% and 
1.04-1.38; 61.10-69.01%, 100%, 97.10-100% and 1.13-1.43; 
66.00-86.61%, 87.87-100%, 90.07-100% and 0.96-1.38, and 
95.14%, 91.05%, 84.91% and 0.94; respectively (Table 2). 
All crosses yielded viable progenies, with no evidence of 
genetic incompatibility being observed among them, ex- 
cept for only sterile Fi-progeny males of which the atro- 
phy of testes and accessory glands were recovered from 
the reciprocal crosses in both directions {An. sinensis x 
An. kleini and An. kleini x An. sinensis) (Figure 1A), while 
normal development occurred in all males from repeated 
backcrosses (Figure IB). On the other hand, all females 
from F x -hybrids and repeated backcrosses yielded normal 
development of ovarian follicles (Figures 1C and D). Re- 
garding repeated backcross groups, the experiments of 
{An. sinensis x An. kleini) x An. sinensis were carried out 
from F!_ 10 , whereas those of {An. kleini x An. sinensis) x 
An. kleini were investigated from F^. The reason for this 
was that the repeated backcross of {An. kleini x An. 
sinensis) x An. kleini could be carried out in only the 
fifth generation, which led to a lack of hybrids for fur- 
ther experimentation (repeated twice). Regarding hy- 
brid crosses, the experiment was performed only on the 
F n of [{An. sinensis x An. kleini) x An. sinensis] x [{An. 
sinensis x An. kleini) x An. sinensis]. 

PCR species identification 

For molecular identification, sequences of nuclear ITS2 
and mitochondrial COI genes were determined for both 
the parent mosquitoes and their hybrid progenies. Se- 
quences of An. sinensis parents (F 0 ) were well matched 
in both genes with those in the GenBank nucleotide se- 
quence database. The four specimens of an An. kleini 
parent (F 0 ) were also checked, and they all matched An. 
kleini based on ITS2 sequences. However, their COI 
gene, An. kleini matched in both species, with 1 being 
An. kleini and 3 An. sinensis (Table 1). 

The F^hybrid progenies have mixed sequences in their 
nuclear ITS2 gene. They have heterogeneous ITS2 se- 
quences because they receive ribosomal RNA genes, in- 
cluding ITS2 from both parents. As for a maternal, 
mitochondrial COI gene, the hybrids followed the trait of 
the mother. All the progenies of An. sinensis x An. kleini 
matched An. sinensis based on COI sequences. However, 
progenies of An. kleini x An. sinensis matched both An. 
kleini and An. sinensis based on COI sequences, in which 
two progenies matched An. sinensis and one An. kleini 
(Table 1). 

Two backcrossing groups were checked, with one be- 
ing a progeny of {An. sinensis x An. kleini) x An. sinensis 
and the other a progeny of {An. kleini x An. sinensis) x 
An. kleini. In the {An. sinensis x An. kleini) x An. 
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Table 2 Hybridization experiments between isolines of An. sinensis and An. kleini 



Crosses (Female x Male) 


Total eggs 
(number)* 


Embryonation 
rate + 


Hatched 

n (%) 


Pupation 

n (%) 


Emergence 

n (%) 


Total emergence n (%) 
Female Male 


Parental crosses 
















An. sinensis x An. sinensis 


538 (258, 280) 


96 


498 (92.57) 


463 (92.97) 


449 (96.98) 


211 (46.99) 


238 (53.01) 


An. kleini x An. kleini 


469 (263, 206) 


89 


413 (88.06) 


322 (77.96) 


303 (94.10) 


136 (44.88) 


167 (55.12) 


Reciprocal crosses 
















An. sinensis x An. kleini 


529 (288, 241) 


86 


450 (85.07) 


441 (98.00) 


436 (98.87) 


253 (58.03) 


183 (41.97) ++ 


An. kleini x An. sinensis 


422 (239, 183) 


83 


350 (82.94) 


350 (100.00) 


350 (100.00) 


1 78 (50.86) 


172 (49.1 4) ++ 


Back crosses 
















(An. sinensis x An. kleini)F ] x An. sinensis 


401 (218, 183) 


63 


245 (61.10) 


245 (100.00) 


238 (97.10) 


140 (58.82) 


98 (41.18) 


{An. sinensis x An. /c/e/'n/') F 2 x An. sinensis 


397 (207, 190) 


88 


332 (83.67) 


309 (93.07) 


309 (100.00) 


173 (55.96) 


136 (44.04) 


{An. sinensis x An. /c/e/'n/') F 3 x An. sinensis 


386 (211, 175) 


74 


282 (73.05) 


265 (93.97) 


265 (100.00) 


143 (53.96) 


122 (46.04) 


(An. sinensis x An. /c/e/'n/') F 4 x An. sinensis 


413 (240, 173) 


79 


319 (77.24) 


284 (89.02) 


259 (91.20) 


137 (52.90) 


122 (47.10) 


(An. sinensis x An. /c/e/'n/') F 5 x An. sinensis 


412 (232, 180) 


69 


282 (68.45) 


282 (100.00) 


254 (90.07) 


139 (54.72) 


115 (45.28) 


(An. sinensis x An. /c/e/'n/') F 6 x An. sinensis 


409 (226, 183) 


86 


318 (77.75) 


308 (96.86) 


295 (95.78) 


171 (57.97) 


124 (42.03) 


(An. sinensis x An. kleini)F 7 x An. sinensis 


365 (147, 218) 


87 


307 (84.11) 


301 (98.05) 


295 (98.01) 


157 (53.22) 


138 (46.78) 


(An. sinensis x An. /c/e/'n/') F 8 x An. sinensis 


355 (154, 201) 


84 


284 (80.00) 


281 (98.94) 


267 (95.02) 


139 (52.06) 


128 (47.94) 


(An. sinensis x An. /c/e/'n/') F 9 x /An. sinensis 


315 (190, 125) 


89 


271 (86.03) 


271 (100.00) 


271 (100.00) 


136 (50.18) 


135 (49.82) 


(An. sinensis x /An. /c/e/'n/')F 10 x /An. sinensis 


336 (162, 174) 


94 


291 (86.61) 


266 (91.40) 


257 (96.62) 


126 (49.03) 


131 (50.97) 


(An. kleini x /An. sinensis)^] x /An. /c/e/'n/' 


497 (239, 258) 


74 


343 (69.01) 


343 (100.00) 


343 (100.00) 


182 (53.06) 


161 (46.94) 


(/\n. /c/e/'n/' x /An. sinensis)F 2 x /An. /c/e/'n/' 


427 (226, 201) 


77 


305 (71.43) 


268 (87.87) 


268 (100.00) 


147 (54.85) 


121 (45.15) 


(/\n. /c/e/'n/' x /An. sinensis)F 3 x /An. /c/e/'n/' 


421 (209, 212) 


70 


286 (67.93) 


257 (89.86) 


257 (100.00) 


126 (49.03) 


131 (50.97) 


(/\n. /c/e/'n/' x /An. sinensis)F 4 x /An. /c/e/'n/' 


458 (211, 247) 


79 


328 (66.00) 


302 (92.07) 


287 (95.03) 


155 (54.01) 


132 (45.99) 


Fn hybrid crosses 
















[(/\n. sinensis x /\n. kleini)F u x /An. sinensis] 


329 (171, 158) 


96 


313 (95.14) 


285 (91.05) 


242 (84.91) 


1 1 7 (48.35) 


125 (51.65) 


x [(/\n. sinensis x /An. kleini)F u x /An. sinensis] 

















*Two selective egg-batches of inseminated females from each cross; + Dissection from 100 eggs; n = number. 
++ Sterile male hybrids with atrophy testes and accessory glands. 



sinensis back crossing group, all the progenies were well 
matched An. sinensis NCBI sequences in the ITS2 se- 
quence, except for 3 of F 2 -hybrid progenies (mixed). Also, 
all sequences of the COI gene showed An. sinensis. In the 
{An. kleini x An. sinensis) x An. kleini backcrossing group, 
all sequences of the ITS2 showed An. kleini, but sequences 
of the COI gene matched An. kleini only in 3 of F 2 -hybrid 
progenies, and the remaining F 3 _ 5 progenies matched An. 
sinensis (Table 1). 

Phylogenetic analysis 

The neighbor-joining (NJ) tree was constructed based on 
COI sequences in order to determine sequence divergence 
among the species examined (Table 1, Figure 2). The NJ 
tree showed concordant results with the PCR assay in all 
samples, which were divided into two major clades with 
67-74% bootstrap support. Clade I consisted of 55 samples 
of An. sinensis. It is interesting to note that among the 55 
samples, 12 (KFO-2, KFO-3, KFO-4, KSKF3-1, KSKF3-2, 
KSKF3-3, KSKF4-1, KSKF4-2, KSKF4-3, KSKF5-1, KSKF5- 



2 and KSKF5-3) and 8 (SKF1-1, SKF1-2, SKF1-3, KSF1-1, 
KSF1-3, SKSF2-1, SKSF2-2 and SKSF2-3), were identified 
as An. kleini and mixed sequences of both species, re- 
spectively, based on ITS2 sequences. Furthermore, these 
samples were placed within the same clade as the pub- 
lished sequence of An. sinensis (mean genetic distances = 
0.003). Clade II comprised 5 samples of An. kleini derived 
from parental (KFO-1), reciprocal (KSF1-2) and back- 
crosses (KSKF2-1, KSKF2-2 and KSKF2-3) based on COI 
sequences (Table 1). The mean genetic distance between 
these two species was 0.023. 

Pyrosequencing analysis 

The 24-nucleotide target region of the COI gene, includ- 
ing positions 237-260 of An. sinensis and An. kleini, was 
useful in classifying the sequenced species into 2 groups, 
as shown in Table 3 and Figures 3 and 4. The first group 
consisted of Parental: An. sinensis F 0 (Figure 4A), F 2 : An. 
sinensis x An. kleini - > hybrid F 2 (Figure 4B), F 5 : hybrid 
F 4 x An. sinensis - > hybrid F 5 (Figure 4C), F 10 : hybrid F 9 
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Figure 1 Reproductive system of adult males and females. (A) F r hybrid adult male of An. sinensis x An. kleini showing atrophy of testes and 
accessory glands. (B) Backcross adult male of {An. kleini x An. sinensis)^] x An. kleini showing normal development of testes and accessory glands. 
(C) F r hybrid adult female of An. sinensis x An. kleini showing normal development of ovarian follicles. (D) F r hybrid adult female of An. kleini x An. 
sinensis showing normal development of ovarian follicles. 



x An. sinensis - > hybrid F 10 (Figure 4D), and F x : An. 
kleini x An. sinensis - > hybrid F 2 (Figure 4F) and an- 
other group comprised Parental: An. kleini F 0 (Figure 4E) 
and F 5 : hybrid F 4 x An. kleini - > hybrid F 5 (Figure 4G). 
Both groups differed from each other in three nucleotide 
positions (T237C, A243G and C253T) (Table 3). The 
positive -control plasmids showed similar results to the 
seven samples in the pyrogram, whereas, a negative con- 
trol did not provide the pyrogram result. 

Discussion 

Introgression or introgressive hybridization is the move- 
ment of a gene (gene flow) between species through 
hybridization by repeated backcrossing of an interspecific 
hybrid with one of its parent species. It can have important 
effects on dynamics of the hybrid zone, speciation and 
adaptive radiation [41]. The variation of mitochondrial 
DNA is important and used widely for indirect studies of 
gene flow [42]. Thelwell et al. [43] reported evidence of 
mitochondrial (ND5) introgression between An. bwambae 
and An. gambiae. Consequently, extensive investigations of 
introgression between An. gambiae and An. arabiensis, and 
An. bwambae and An. gambiae have been documented sys- 
tematically and extensively during the past decade [44-46]. 
Additionally, Walton et al. [42] demonstrated that mito- 
chondrial DNA introgressed from An. baimaii (= dims D) 
into An. dims (= dims A), and Morgan et al. [47] reported 
evidence to support the mitochondrial introgression 



between An. baimaii and An. dims by the high levels of bi- 
directional mitochondrial gene flow detected between these 
2 species. A large number of anopheline species were re- 
ported to be capable of interspecific hybridizations under 
laboratory conditions [48]. However, few species have suc- 
ceeded in natural hybridization, for example, between 
An. gambiae and An. arabiensis [23], An. bwambae and 
An. gambiae [43], An. scanloni (= dims C) and An. baimaii 
(= dims D) [24], and An. minimus and An. harrisoni 
[28,29]. Recent discovery of the natural hybrid between An. 
sinensis and An. kleini [27], and successful establishment of 
iso-female line colonies of these 2 anopheline species has 
urged this study to form a reliable systematic procedure to 
confirm this natural event and/ or perform an introgressive 
study. The results of sterile F 2 -hybrid adult males with atro- 
phy of accessory glands and testes obtained from the recip- 
rocal crosses between these 2 anopheline species are in 
keeping with "Haldane Rule", which states that in interspe- 
cific crosses, the heterogametic sex (X, Y) will show sterility 
or viability problems before the homogametic sex [49]. 
Thus, the results of this study agree with those of Davidson 
[50], who reported that female and male hybrids obtained 
from An. gambiae and An. arabiensis were fertile and ster- 
ile, respectively. Furthermore, the results of this study are 
in accordance with crossing studies in the laboratory by 
Baimai et al. [51]. They demonstrated that F 2 hybrids ob- 
tained from the cross of female An. dims with male An. 
scanloni were fertile and viable, with an exception of sterile 
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KFO-2 
Anopheles sinensis 




HF11-2 
SKSF9-3 

SKSF2-3 



Figure 2 Neighbor-joining tree of /An. sinensis and /An. /c/e/n/ based on COI sequences. Numbers on branches are bootstrap values (%) of NJ 
analysis. Only greater than 50% bootstrap values are shown. Bars represent 0.005 substitutions per site. Detailed code of samples is shown in Table 1. 



Table 3 Sample details and nucleotide positions used for discriminating between An. sinensis and An. kleini based on 
COI sequences 



Code no. 


Samples 
(Female x Male) 


237 


Nucleotide at position 
243 


253 


A1* 


Parental: An. sinensis F 0 


T 


A 


C 


A2* 


Fv An. sinensis x An. kleini - > hybrid Ft 


T 


A 


C 


A3* 


F 5 : hybrid F 4 x An. sinensis - > hybrid F 5 


T 


A 


c 


A4* 


F 10 : hybrid F 9 x An. sinensis - > hybrid F 10 


T 


A 


c 


A5 + 


Parental: An. kleini F 0 


C 


G 


T 


A6* 


Fv An. kleini x An. sinensis - > hybrid Ft 


T 


A 


C 


A7 + 


F 5 : hybrid F 4 x An. kleini - > hybrid F 5 


C 


G 


T 



*group 1, COI sequence matched with An. sinensis. 
+ group 2, COI sequence matched with An. kleini. 
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A1 
A2 
A3 
A4 
A6 
A5 
A7 

consensus 



GflGCCCCTGflTflTflGCTTTTCCT 
GflGCCCCTGflTflTflGCTTTTCCT 



GHGCCCCTGflTfiTflGCTTTTCCT :GHRTRRHTfiR 
GflGCCCCTGflTflTflGCTTTTCCT :GflflTflflHTftfl 
GflGCCCCTGflTflTflGCTTTTCCT :GARTRAATAR 
GRGCTCCTGRTflTRGCTTTTCCT :GRRTRRRTRR 
GflGCTCCTGflTflTflGCTTTTCCT IGRRTRRRTRR 
GflGCcCCTGRTflTRGCTTTTCCT|GRRTRRRTRR 

Forward primer 



GRRTRRRTRR 
GRRTRRRTRR 



231 

flTRflGTTTTTGflflTflTTflClTCCTTCflTTfiftCflTTflCTflflTTTC 
ATAAGTTTTTGflflTRTTRC JCCTTCflTTRflCflTTflCTfiflTTTC 
UTARGTTTTTGflflTATTflC TCCTTCflTTBflCflTTflCTflflTTTC 
ATflHGTTTTTGflflTATTflC TCCTTCflTTAflCATTflCTAATTTC 
flTflflGTTTTTGflflTflTTflC ITCCTTCflTTAftCflTTACTAATTTC 
flTflflGTTTTTGflflTflTTflC CCCTTCGTTAflCflTTATTflATTTC 
flTflflGTTTTTGflflTflTTflC ECCTTCSTTAflCATTATTAflTTTC 



RTRRGTTTTTGRfiTRTTf 

Sequencing primer 



AC JCI 



CTTC^TTftflCflTTfl^TflflTTTC 

Target region 



AGAAGT ATAGT AGAAAATGGAGCAGGA AC AGGGTGAAC T G 1 TTflTCCCCCflCTTTCflTCTCG 
RGRRGTflTflGTflGRflRRTGGflGCflGGflRCRGGGTGRflCTGl rTRTCCCCCftCTTTCflTCTGG 
nGAnGTflTnGTflGflflnflTGGflGCflGGflncnGGGTGflRCTGl TTflTCCCCCftCTTTCflTCTGG 

AGAAGTATAGTAGAAflflTGGflGCfl6GflflCAGGGTGAACTGl TTflTCCCCCflCTTTCflTCTCG 
flGflflGTflTflGTAGflflflflTGGAGCRGGflflCflGGflTGflRCTGl ITRTCCCCCTCTTTCflTCTGG 

RGAAGTATAGTAGAAAATGGAGCAGGAACAGG 8 TGAACTG1 TTflTCCCCCaCTTTCflTCTCG 

Reverse primer 



Figure 3 Alignment of the cytochrome c oxidase subunit I (COI) gene derived from seven samples. Position of the forward primer 
(Anop_COI_F) and biotinylated reverse primer (Anop_COI_Rb) for template amplification are shown in the black boxes, while the sequencing 
primer (Anop_COI_S) and target region are shown in the green and red boxes, respectively. Asterisks indicate position of the target region used 
for species level identification. 



males. It was interesting to note that the repeated back- 
cross progenies in both directions, which resulted in 
obtaining An. sinensis from hybrids of F 2 -n progenies 
[{An. sinensis x An. kleini) x An. sinensis] and F 3 - 5 pro- 
genies [{An. kleini xAn. sinensis) xAn. kleini], indicated 
the presence of introgressive hybridization between An. 
sinensis and An. kleini. 

Regarding PCR identification of parental specimens, 
the exact species of An. sinensis and An. kleini were 
used in this study based on both the nuclear ITS2 and 
mitochondrial COI genes. Interestingly, 4 iso-female 
lines of An. kleini showed the correct gene trait in the 
ITS2, but their COI sequences matched both species (i.e., 
An. kleini: 1 iso-female line, and An. sinensis: 3 iso-female 
lines). Therefore, the authors assumed that the COI se- 
quences of An. kleini had been replaced by those of An. 
sinensis. Subsequently, reciprocal and repeated back- 
crosses were performed to clarify our hypothesis. Most 
COI sequence results of the progenies obtained from re- 
ciprocal and repeated backcrosses revealed that the 
mitochondrial COI gene introgressed from An. sinensis 
into An. kleini. This event resulted from introgression 
that occurred between these two species via consecutive 
repeated backcrosses. Thus, the sample identification of 
these 2 wild-caught species should be careful only when 
the COI barcoding region has been applied. Also, re- 
sults from phylogenetic analysis confirmed the existence 
of an introgression phenomenon between them. Fur- 
thermore, the results from this study are in agreement 
with those of Petit and Excoffier [52]. They suggested 
that in species with male-biased dispersal (heterogam- 
etic sex), mtDNA markers should introgress more read- 
ily than nuclear ones. 

Pyrosequencing is a unique sequencing method that was 
developed as an alternative to classical DNA sequencing 
for short- to medium-read applications. It is an accurate, 
simple and flexible bioluminometric method, which does 
not need labeled nucleotides or gel electrophoresis [53]. 
To date, this technology has been used successfully for 
high throughput identification of bacteria [54,55], virus 
[56,57], protozoan parasites [53,58-60] and helminthes 



[40]. This study applied the pyrosequencing technology 
for reliable identification from seven samples of An. sinen- 
sis and An. kleini. The results were consistent with those 
of dideoxy sequencing and phylogentic analysis in six sam- 
ples, except for the one (sample code no. A7) matched 
with An. kleini COI sequences. This sample possibly had a 
different mitochondrial haplotype from other hybrid pro- 
genies. More recently, the next-generation sequencing 
provided a good explanation of interspecific gene flow be- 
tween An. gambiae and An. arabiensis [61] and An. gam- 
teMandS [62,63]. 

The effectiveness of a vector control strategy, and genet- 
ically modified strains of mosquitoes in a population that 
is unable to transmit malarial parasites, relies upon the 
gene flow within species and introgression [64]. Remark- 
ably, Rheindt and Edwards [65] mentioned concern that in 
the long-term introgression with newcomer species may 
lead to a loss of genetic integrity in native species. Detec- 
tion of natural hybridization is rare but meaningful in 
terms of horizontal transfer of advantageous genes, such 
as those in malaria susceptibility, particularly when species 
that allow hybridization are susceptible or refractory to 
malarial parasites [9,66,67]. This event is also involved in 
insecticide resistant genes, e.g., the kdr gene [68,69] and 
ace-1 gene [70]. In addition, Morgan et al. [71] stated that 
the absence or presence of gene flow between populations 
and species has an impact on the dynamics of malaria 
transmission as well as construction of effective strategies 
for controlling malaria vectors. Our studies presented 
the introgressive events through consecutive repeated 
backcrosses under laboratory conditions, in which the 
mtDNA gene could be moved from one species to an- 
other. However, natural movement of the refractory 
gene to vivax malaria between sympatric populations of 
a low potential vector {An. sinensis) and a high potential 
vector {An. kleini) needed intensive and systematic 
clarification. This new body of knowledge is anticipated 
to elucidate the promising strategies for replacing popu- 
lations of high potential vectors with that of low poten- 
tial vectors by using genetic manipulation through the 
gene introgression mechanism. 
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Dlsp: A TC G AT C GAT CGA T C GA TCG AT C GAT CG A TC G 



Dkp: A TC G AT C GAT CGA T C GA TCG AT C GAT CG A TC G 



125 
120 
115 
110 



XL 



125 
120 
115 
110 



E S A TC GA TC G AT CG AT C GA TC GA T CG AT CG A TC G 
I 5 10 15 20 25 30 

Scq: TCC TTC ATT AA C ATT ACT AATTTC 

1 llllnlilll. 

Dlsp: A TC G AT C GAT CGA T C GA TCG AT C GAT CG A TC G 



E S A TC GA TC G AT CG AT CG A TCG A TC GAT C GA TC G 
■ 5 BO 15 20 25 30 

Seq: TCC TTC ATT AA C ATT ACT AATTTC 




Dlsp: A TC G AT CGA T CGA T C GA TCG AT C GAT CG 



120 
115 
110 
105 





















1,1; 




1 I 


i..y 1 





130 
120 
110 



JLL 



JJJLi 



Ik 



E S A TC GA TC G AT CG AT CG A TCG A TCG AT C GA TC G 
I 5 10 15 20 25 30 

Scq: TCC TTC ATT AA C ATT ACT AATTTC 



E S A TC GA TC G AT CG AT CG AT CG AT C GAT CG 
I 5 110 15 20 I 25 

Scq: CCC TTC G TT AA C ATT A TT AATTTC 




1 1 31 I Lll 1 1 I 

Dlsp: A TC G AT C GAT CGA T C GA TCG AT C GAT CG A TC G 



Dlsp: A 




TCGATCGATCGATC GAT CG A TC G 



125 
120 
115 
110 



125 
120 
115 
110 




E S A TC GA TC GAT CG AT C GA TC GA T CG AT CG A TCG 
I 5 ■10 15 20 25 30 

Seq: TCC TTC ATT AA C ATT ACT AATTTC 

G 



ES A TC GA TCGATCG ATCG ATCG A TCGATC GATCG 

I 5 10 15 20 25 30 

Seq: TCC TTC ATT AA C ATT A C T AATTTC 




I II IE III 

Dlsp: A TC G AT CGA T CGA T C GATCG AT C GAT CG 



111 



125 
115 
105 



JLU. 



ES ATCGATCGATCG ATCG ATCG ATCGATCG 

5 10 15 20 25 

Seq: CCC TTC G TT AA C ATT ATT AATTTC 

Figure 4 Pyrograms showing sequence analysis (SQA) of 24-base fragments of the COI gene. Group 1; (A) Parental: An. sinensis F 0 , (B) 
An. sinensis x An. kleini - > hybrid F 1; (C) F 5 : hybrid F 4 x An. sinensis - > hybrid F 5; (D) F 10 : hybrid F 9 x An. sinensis - > hybrid F 10 , and (F) Fy.An. kleini x An. 
sinensis - > hybrid F v Group 2; (E) Parental: An. kleini F 0 and (G) F 5 : hybrids F 4 x An. kleini - > hybrid F 5 using pyrosequencing. Theoretical pyrogram 
patterns (top of each panel) and representative raw data (bottom of each panel) of control DNA extracted from each An. sinensis and An. kleini by 
pyrosequencing are shown. Pyrosequencing was performed by addition of enzyme (E), substrate (S), and four different nucleotides. The letters under 
the black bars show the dispensation (Disp:) order. The actual sequence detected by pyrosequencing is indicated below the panels after "Seq". The 
Y-axis represents the level of fluorescence emitted by incorporating a nucleotide base, and the X-axis represents the total number of bases added at 
that point in time; A, T, C, G nucleotide bases. The light gray areas show the pyrogram for identifying each of the two groups of An. sinensis and 
An. kleini. 
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Conclusions 

A single gravid, natural hybrid female between high {An. 
kleini) and low {An. sinensis) potential vectors of P. vivax 
was discovered in Paju City, Republic of Korea (ROK). The 
discovery of natural hybrids between these two anopheline 
species has led to systematic investigations of various as- 
pects that clarify this event. Hybridization experiments 
used iso-female line colonies of these anophelines together 
with DNA analysis of ribosomal DNA [second internal 
transcribed spacer (ITS2)] and mitochondrial DNA [cyto- 
chrome c oxidase subunit I (COI)] of the parental colonies, 
Fi -hybrids and repeated backcross progenies were per- 
formed intensively by using a PCR-based assay and py- 
rosequencing technology. The results revealed that 
introgression of the COI gene between An. sinensis and 
An. kleini was involved in this phenomenon. The pure 
An. sinensis obtained from hybrids of repeated back- 
cross progenies in both directions, i.e., F 2 _n progeny 
[{An. sinensis x An. kleini) x An. sinensis] and F 3 _ 5 pro- 
geny [{An. kleini x An. sinensis) x An. kleini] provided 
obvious supportive evidence. The results emphasize a 
promising way to replace the population of a high po- 
tential vector {An. kleini) with that of a low potential 
vector {An. sinensis) through the mechanism of gene 
introgression. 
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